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Abstract 

We compute the CP violating decay asymmetries relevant for baryogenesis sce- 
narios involving the out of equilibrium decays of heavy particles, including the 
finite temperature effects arising from the background of light thermal particles 
which are present during the decay epoch. Thermal effects can modify the size 
of CP violation by a sizeable fraction in the decay of scalar particles, but we find 
interesting cancellations in the thermal corrections affecting the asymmetries in the 
decays of fermions, as well as in the decay of scalars in supersymmetric theories. 
We also estimate the effects which arise from the motion of the decaying particles 
with respect to the background plasma. 



1 Introduction 



A classic, and still very attractive, scenario for the generation of the baryon asymmetry 
of the Universe, is based on the fact that very massive particles fall out of equilibrium 
as the temperature of the Universe drops below their mass and their equilibrium density 
becomes Boltzmann suppressed. If they decay at that epoch through baryon (B) violating 
channels which also violate CP, a net baryon asymmetry will result. 

Usually, in these scenarios the CP violation results from a one-loop decay diagram, 
whose interference with the tree level process allows the phases of the complex coupling 
constants to show up in the decay asymmetries. 

In addition to the complex couplings, to have a non-zero partial decay rate asymmetry 
requires the loop integral to develop an absorptive part. This actually happens whenever 
the loop diagram can be cut in such a way that the particles in the cut lines can be 
produced on-shell, and hence these particles need to be lighter than the decaying one. 
In the proposed baryogenesis scenarios of this kind, the light particles in the loop are 
just standard quarks, leptons or Higgs bosons, so that the appearance of a non-vanishing 
absorptive part is guaranteed. 

At the high temperatures at which the heavy particles decay, the light standard 
particles are in equilibrium with the hot plasma present, and hence a question arises on 
whether the existence of the background particles has any effect in the evaluation of the 
CP violating asymmetries. Indeed, some time ago Takahashi [|^] showed that the thermal 
effects could modify the predictions for baryogenesis in SU{5) models by up to ~ 40% 
with respect to the T = results, and hence these effects may need to be taken into 
account in the proper computation of the resulting baryon asymmetry in specific models. 
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The main effect of the background can be taken into account by employing finite 
temperature propagators in the computation of the loop. In tliis way it is possible to 
consider simultaneously both the 'direct' propagation of a particle between two vertices 
in the loop and the absorption by the medium of a particle from the first vertex combined 
with the emission of another one towards the second. These two alternatives are actually 
indistinguishable in a thermal bath. 

Another implication of the finite density background is the modification of the final 
state phase space density distributions, which take into account the stimulation of the 
decays into bosons and the Pauli blocking of the decays into fermions, and this may 
eventually also affect the rates. 

2 SU{5) triplet decays 

In order to discuss these issues, let us start by reanalysing the scalar decays into two 
fermions, which is relevant in the case of heavy Higgs boson triplet decays in SU (5). We 
will generalise the computation in ref. lU by including also the CP violating diagrams 
arising from mixing among different heavy states P, ^ ^ . Notice that in order to have 
CP violation at one loop in the SU{5) GUT it is necessary to have more than one heavy 
state, and also the existence of more than one heavy state is natural in other models, 
such as the leptogenesis scenarios to be discussed later, and hence the contribution from 
the mixing among the heavy states is generally important. For simplicity we will consider 
the case in which the masses of the heavy states are significantly splitted, since the study 
of the CP violation in the near degenerate situation present additional complications 

One generally expects that only the decay of the lightest of the heavy states, Ti, will 
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be the one leading to a net B, since any asymmetry produced at earlier times through 
the decay of heavier states would be erased by B violating processes which could still be 
in equilibrium^ Hence, one has to compute the asymmetry resulting from the Ti decay, 
and in the case we will consider in which there is a hierarchy among the masses of the 
heavy states (M^ ^ Mi, with /c > 1), it will be natural to assume that the heavier states 
Tfc have already decayed, and hence have a negligible density, at temperatures T < Mi 
when the lighter one is falling out of equilibrium. 

Let us consider the SU{5) lagrangian involving several scalar five-plets $j = (Tj, Hi), 
containing, together with Higgs doublets Hi, the heavy color triplets Tj : 

C = (fr"^ XP) + J^^a (e"^"'^ ^5e) + h.c, (1) 

where the gauge indices are denoted by greek letters. The matter fields are in the decuplet 
and the fiveplet representations as usual, \E' = {q,u'^,e^) and x = {d,l^). Since they are 
vectors in flavour space, the Yukawa couplings and Qi should be thought as 3 x 3 
matrices, but for simplicity the flavour indices are not displayed. 

The CP violating B asymmetry arising from the decay of a Ti and Ti pair is 



_ Y.fBf[V{Ti^Fj)-V{f,^Ff)] 

with Bf the baryon number of the final states Ff = qi"^, ue'^, u'^d, (fq. From the 
interactions in eq. (|l|), at zero temperature we have 





— Im|Tr 


[ghi /./!)} Im{/,(x,)} + ^Im{Tr 




}lm{/,} 




4Tr(/i//) 


+ 3Tr((7l(7i) 





where is the heavy-to-light ratio of squared masses: 

Mf 



(3) 



^However, it is possible to imagine the situation in which a heavier triplet has the smallest decay 
rate, if the effect of the larger phase space is compensated by the smallness of the relevant couplings. 
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Notice that the asymmetry is null if there is a single triplet field. In this case, the result (^) 
depends on the hermitian matrices glgi and fifl; therefore their trace has no imaginary 
part, and the trace of the product is equal to half the trace of the anticommutator, again 
an hermitian matrix. The terms in the denominator correspond to the tree level decay 
rate; those in the numerator to the interference between the tree level and the absorbptive 
part of the loop amplitudes^. In particular. If arises from the loop integral in the 'vertex' 
contribution, in which the is exchanged in the t-channel (see fig. la), and Is comes 
from the loop integral in the 'wave' contribution, in which the Tk is exchanged in the 
s-channel (fig. lb), i.e. 

u{pi)v{p2)It{xk) = u{pi) J J^:^S{p - q, 0)S{-q, 0)D{q - p2, Mk)v{p2), (5) 

/ 7SlTr(5(-g,0)5(p-g,0)), (6) 



^ 2 J {27iy 



with S{p,m) {D{j),m)) being the propagator of a fermion (scalar) with mass m and 
momentum p. 



From these expressions one gets the following absorptive parts: 

Im{/t(x)} '"^ - /- 1 



167r 
and 



1 - a; In ( 1 + 

X 



(7) 



Im{/J ^ 



IGvr 

We consider now the finite temperature effects on the CP violating asymmetry e. The 
main effect comes from using instead of the usual T = propagators, the finite T ones. 
We choose to work in the real time formalism (RTF) of thermal field theory |]I0|. The 



^The apparent difference with respect to refs. |i[ in the denominator is just due to a different 
normahzation of the field 5", and hence of the Yukawa couplings fi. However, our result for the numerator 
differs from ref . in the sign of the first term (in agreement with |^ ) and a factor 2 in the second 
one. 
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Green's functions computed in this formalism are directly the time ordered ones, unlike 
in the imaginary time formalism where different analytical extensions to real momenta 
lead to different Green's functions (retarded, advanced, etc.) [O. The RTF involves 



the introduction of a ghost field dual to each physical field, which leads to a doubling 
of the degrees of freedom. The thermal propagator has then a 2 x 2 matrix structure: 
the (11) component refers to the physical field, the (22) component to the corresponding 
ghost field and the off-diagonal (12) and (21) components mixing them. However, since 
we are working only at one loop and the external legs are physical, i.e. type-1 fields, 
we just need the (11) component of the propagators, which are, for fermions and bosons 
respectively. 



Snip, m) = {■^+ m) 



— - 2TmF{p ■ u)5{j)^ - 

P'^ — + «0+ 



(9) 



Dn{p,m) 



+ 27mB{p ■ u)5{p'^ — 



p2 _ ^2 _|_ 

with u the 4-velocity of the medium [u = (1, 0, 0, 0) in the medium rest frame), and 



(10) 



1 



exp {\x\ / T) ± 1 
We drop the (11) sub index in the propagators from now on. 



(11) 



As discussed previously, we neglect the background density of the heavy triplets Q 
Tfc (assuming ^ Mi), and as a first approximation we assume that the decaying 
particle Ti is at rest (particle motion effects will be considered in section 4). Using the 
well known property, 

= (12) 



^Notice that the exchange of Ti in the one-loop diagrams does not contribute to e, so that only A; > 1 
are relevant. 
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we get for the absorptive part of the vertex loop integral 



Im{/,^(a;,)} = Im{/,(x,)} [l - 2n, + 2^] -^T -^n, (^) , (13) 



nF,B = nF,B ( ^ ) ■ (14) 



where 

- _ /- 

This is similar to the result in ref. except for the relative sign between the T = 
part and the finite temperature corrections, implying that the temperature effects tend 
to reduce (instead of enhancing) the CP violation. The integral term in the r.h.s. of 
eq. (p!3D arises from the absorption of particles from the background which are energetic 
enough so as to put the intermediate state Tk on-shell and hence make it contribute 
to the absorptive part. This term becomes then extremely small in the case in which 
the heavy masses have a significant hierarchy, i.e. for M^, ^ Mi, since the amount of 
background particles which are energetic enough is Boltzmann suppressed. 

For the finite temperature 'wave' contribution we obtain 

Im{/J} = Im{J,} [l - 2nF + 2n|] , (15) 

so that the temperature dependence is similar to the one in the leading term of the vertex 
part. To give a quantitative idea of the effect, we notice that if the temperature is taken 
to be 1, 1/3 or 1/10 of the lightest triplet mass. Mi, the overall factor in ( plSf ) including 



the temperature effects is 0.53, 0.70 or 0.99 respectively. The physical interpretation of 
this effect is simple: due to the thermal background, the two light fermions exchanged 
in the loop are subject to a Pauli blocking, and this leads to a reduction of the amount 
oi CP violation. 



^We have checked that the same resuh can be obtained by applying the Cutkosky cutting rules at 
finite temperature in the real time formalism pM- 
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In this scenario there is no effect on e resulting from the final state blocking, since this 
just leads to overall factors (1 — ni?)^ multiplying the rates, and hence these factors cancel 
in the ratio in eq. (^. Other thermal effects, such as thermal masses or wave function 
renormalization, are higher order in the coupling constants and hence we neglect them. 

3 Leptogenesis scenarios 

The SU (5) model discussed in the previous section, in spite of being the prototype for the 
'out of equilibrium decay' scenarios of baryogenesis, has the drawback that it generates 
no net B — L asymmetry (a characteristic of SU{b)), and hence the B generation is 
vulnerable to the anomalous B violating processes of the Standard Model |]13[ (which 
only leave B — L unaffected) , with the consequence that all asymmetries generated within 
this model will be eventually erased. 

A very interesting way out to this problem [0, |T5[ is based on the generation of a 
lepton (L) asymmetry at early times, by the out of equilibrium decay of heavy isosinglet 
neutrinos (the usual ones appearing in see-saw models for neutrino masses and naturally 
present in GUT models such as 50(10)). 

In this section we discuss temperature effects in this kind of models, first under the 
minimal assumption that the standard model spectrum is enlarged to include heavy right- 
handed neutrinos, and then, considering the supersymmetric extension of the model. 

3.1 Non- supersymmetric case 

The interactions of the heavy neutrinos A^j, in the basis in which their mass matrix is 
diagonal, are given by the following Lagrangian 

C = -Xa^e^pN,PdaH^ + h.c. (16) 
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where la = (i^a, la) H = {H^, H^) are the lepton and Higgs Standard Model 
doublets (a = e, fi,T, i = 1, 2, 3, and e^/? = — C/sq, with eu = +1). 

The complete T = CP violating L asymmetry was computed for this model in 
ref. 0, resulting in 



k>l 



Im{ Jt(xfc)} + 2^ — -Im{ J,} 



Xk-l 



(17) 



where Xk is defined analogously to eq. (HI) and 



Im|(AtA) 



(AtA)n ■ ' ' 

Notice that, in full analogy with the SU{5) case, a single right handed neutrino would 

be unable to generate any asymmetry. 



The loop integrals are given in this case by 



uipi)PRuip)Jtixk) = u{pr)PR / j^^S{q, 0)5(g - p2, Mu)D{p - q, O)PiM(p) 



(2vr) 



and 



d^q 



u{pi)u{p)MiJs = -iu{pi) I j—^S{q,0)D{p- q,0)u{p), 



so that the absorptive parts result 



Im{Jt(x)} 



IGtt 



1 - (1 + a;)ln 1 + - 

X 



and 



Im{Jj 



327r 



(19) 



(20) 



(21) 



(22) 



The computation of the finite temperature contribution to these quantities is similar 
to the one in the previous section, and leads to 

Imj Jf (xfc)! = lm{Jt{xk)} [1 - np + Ub - 2nFnB] 



y/x^ f°° du 



167r 



Up 



/Miu 

Im| jjj = Im{ J^} [1 - Up + ub - 2nFnB\ , 



du 
Jxk u + 1 



[U - Xk) + Ub 



/Miu 



{Xk + 1] 



(23) 



where for instance in this last expression the different contributions in the r.h.s. clearly 
separate into the pieces coming from the T = propagators, the one from the thermal 
correction to the fermion {£) propagator, the one from the thermal piece of the boson 
{H) propagator and the product of these two corrections. However, it is easy to check 
that the Bose and Fermi distributions satisfy 

nsiE) - npiE) = 2nB{E)nF{E), (24) 

so that the main temperature correction cancels out (only the small integral term in 
Im{ Jj} survives). This is due to the opposite effects resulting from the Pauli blocking of 
the loop fermion line and the stimulation of the bosonic loop line. On the other hand, 
here again the final state statistical factors {1 — np){l + ub) cancel in the expression for 
e, and as a result no significant temperature dependent effect is found. 

3.2 Supersymmetric case 

This scenario has also received considerable attention within a supersymmetric framework 
p!6| , 0, in particular because the scalar partner of the heavy neutrino A'^i is a good 
candidate for being the inflaton field, in which case the L asymmetry could be produced 



during the process of reheating of the Universe as A^^i decays fl^. In this case, as we will 
now show, another interesting cancellation is found in the asymmetry produced by Ni 
decay. 

The scalar neutrino can decay either into two scalars {LH) or into two fermions (ih), 
and the contribution to CP violation in one channel is obtained from the loop involving 
the particles of the other channel P]. The thermal effects will modify the asymmetries 
corresponding to each channel, in a way similar as they did in the case of SU{5). For 
instance, in the LH channel, if we ignore the small integral piece coming from the vertex 
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(equivalent to the last term in the r.h.s. of eq. (0)), the asymmetry will be 

e^{N LH) ~ e^=°(iV ^ LH) [1 - 2^^.(1 - np)] , (25) 

and similarly 

e^(iV fh) ~ e^=°(iV ^ [1 + 2nB(l + n^)] . (26) 

However, due to the effects of the final state phase space factors (1 =F ?^f,b) entering 
into the partial decay rates, the branching ratios of the two different channels will no 
longer be equal (as is the case at T = 0). One has instead that 

BRiN ^ih)= ^^_ - 1 - BRiN LH). (27) 

The total asymmetry produced in the N decay is 

= BR{N ^ fh)e^{N ^ ih) + BR{N LH)e^{N LH), (28) 

with the surprising result that the main corrections arising from thermal effects actually 
cancel out, leading to 

4 - 4^'- (29) 

Notice that there are also new supersymmetric diagrams contributing to the CP 
violating asymmetry in the heavy neutrino decays. However, since the particles in the 
loop as well as the external ones are always one fermion and one boson, both massless, 
the cancellation found in the previous subsection will also occur in the new channels, and 
therefore there are no significant thermal corrections to the zero temperature result of 
ref. 0. 



4 Effects of particle motion 



In the discussion so far we have always considered the decay rate of a particle at rest in 
the thermal bath. However, the decaying particle will in general be moving through the 
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background with non zero velocity (3 = v/c. Since now the Lorentz symmetry is exphcitly 
broken by the plasma, this motion can in principle affect the thermal corrections to the 
decay asymmetry. When the leading thermal corrections cancel, as in the leptogenesis 
scenarios discussed in Section 3, the effects of the motion of the decaying particle will 
provide the main thermal corrections, and hence it is worth to quantify them. To estimate 
the size of these effects, we will consider here the case of the heavy neutrino decay in the 
non-supersymmetric model. The other cases can be analysed similarly. 

It is convenient to compute the decay rate asymmetries in the rest frame of the 
decaying particle, where the medium will be characterised by a non trivial 4-velocity 



u = (7,-7/5), with 7 = l/vl— 7^ as usual. In this system, the effect of the motion 
will reflect in a modification of the equilibrium distributions appearing in the thermal 
propagators, eqs. (pHTOj). The decay rate will also depend on /? through the final state 
phase space factors, which for the case of fermion decays is 

[1 - riF^h ■ u)] [1 + nB{kH ■ u)] = ^_^^p^^^^^^^^^ P;3(cos^), (30) 

where 6 is the angle between f3 and the momentum of the final state lepton, and we have 
introduced 

We see that now there is a privileged direction selected by the plasma spatial velocity 
and therefore the decay process is no more isotropic: the rate depends on the angle of the 
decay product trajectories with respect to this direction. Due to statistics the fermions 
(bosons) are preferentially emitted in the direction parallel (antiparallel) to the plasma 
velocity, which corresponds to the less (more) occupied region of the thermal distribu- 
tion. Anyway, since we are interested in the total decay rate, the angular dependence is 

integrated out and we are left only with a dependence^. 

^The integrated decay rate can depend only on the Lorentz invariants — Aff and p-u — Mij [ p^ . 
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Let us define e/3 as the integrated asymmetry generated in the decay of a heavy 
neutrino moving with velocity (3. The decaying particles will actually have a distribution 
of velocities with occupation numbers n{E), where E = M17. To estimate the overall 
effect of the particle motion we may just approximate this distribution with the Fermi 
Dirac distribution, n{E) = np{E), and compute the average asymmetry 

where A^^i is the particle's volume density and 

i 

,/3Y n(7Mi). (33) 



dm Ml ^2 5 



dl3 271^' 

It is clear that (e), computed for a given temperature T, is just the asymmetry that 
would result if the initial thermal population of heavy neutrinos went out of equilibrium 
and decayed all simultaneously at temperature T. We will use the asymmetry (e) as an 
indicator of the possible effects of the particle motion, although clearly to obtain the 
exact impact of this into the final lepton asymmetry would require to integrate the whole 
Boltzmann equations, a task beyond our scopes. 

To obtain the decay asymmetry at a fixed velocity e/3, it is convenient to separate the 
angular dependences arising from the final state phase space factor Pp{cos9) and the one 
arising from the one loop integrals (via the anisotropic background density in the particle 
rest frame). This last will be included in the factor LpicosO^Xk) defined through 

4 = TT-t^TTTtE^^i r dzLl{z,Xk)Pp{z). (34) 

Here the supraindex c labels the two different contributions to the CP asymmetry in the 
decay, i.e. the "wave" (e*) and the "vertex" (e*) pieces. Clearly at zero velocity we have, 
according to eq. (|lj 



Ll{z,Xk) = ^^Im{Jj, (35) 

Xk i 

Ll{z,Xk) = Im{Jt(xfc)}, (36) 
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which are actually independent of z as expected. 



At /? 7^ we find after direct computation of the interference terms 



p{i)i 



(37) 



where the loop functions /, 



(n) 



Jf3 



dy P(3iy) 



(3^ 



arise after integrating over the angle of the momenta in the loop. Notice that the phase 
space factor appears also in these integrals, arising from the statistical factors in the 
thermal loops. In the limit of zero velocity, we have /d°'' = 1 and Jq^^ = 0. 



Now we can write in a simple form the final expression for ei; namely, putting together 



eqs. (PD, (M) and m) we get: 



'Ni Jo 



dp 



dNi (0) 



dl3 



J a 



f, 



(0) 

/3 



(39) 



It is not possible to integrate analytically this expression, but expanding the loop 
functions to first order in (3"^ we obtain 

2 



1 + 7T 



Mi/2T V 



3^ Vsinh(Mi/2T)y 



(1) 

13 



1 Mi/2r 

3 ^ sinh (Mi/2T) ' 



(40) 
(41) 



Substituting in eq. (pOf), we get 



9^^ ^ Uinh(Mi/2T) 



(42) 



For Mi/T = 1,3,10 we have (e'^)/eg ~ 1.168,1.057,1.000; we see therefore that the 
effect is small compared to the /3 = piece. In the previous evaluation we have used 
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for the average velocity just a simple approximation, which is accurate to the 10% level, 
writing {(3"^) ~ {l3)%oitzi with the average velocity for a Boltzmann distribution of massive 
relativistic particles being (with x = M/T) 

2(1 + x)exp(— x) 



Boltz 



x'^Koix) 



(43) 



Computing numerically eq. (|39|) we find, for the same reference values of Mi/T as before, 
that (e*)/eQ = 1.231, 1.054 and 1.000, so that the approximate result is acceptable. 

Let us now consider the "vertex" contribution. In this case the loop integration is 
more involved and we get the following expression (neglecting the small integral piece 
which is Boltzmann suppressed): 



L^f^{z, Xk] 



IGtt 



4°^ -9p{xk,z) 



where 



1 •l[y + z{l + 2xuW + 4xfc(xfc + 1)(1 - z'') 



(44) 



(45) 



For /5 = the function g actually does not depend on z and is 



go{xk) = (1 + Xfc)ln IH 

Xk 



(46) 



so that we recover exactly eq. (|3q). 

Substituting the result of the loop integration in eq. ( P^ we finally obtain 



kl 
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(0) 

/3 



!-ld'Z 9f3iXk,z)Pf3{z) 



(0) 



(47) 



Again we can evaluate this function analytically only for small /?, by using the expansion 



M /IT I 
9,{xk.z) ^ 9oixk)+Pz-^^^^^^h{xk)+P'[^^^ 



h(Mi/2T) J 



j{xk) - zH{xk) , (48) 
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where 



h{x) = 2{l + x) - (l + 2x)go{x), 

j{x) = (1 + x){l + 2x) -2x{l + x)go{x), 

l{x) = 3{l + x){l + 2x)-{l + 6x + 6x^)go{x). 



(49) 
(50) 
(51) 



So, the average CP violation is given by 



Mi/2T y/ 1 
sinh (Mi/2T)j \gQ{xk) 



- 2xi 



(52) 



We then see that the effect due to the particle motion is again to increase the vertex 
CP asymmetry, as in the case of the wave part. However, the effect depends also on x^, 
i.e. on the hierarchy between the particle masses. In the limit of large Xk we can use the 
expansion {go^Xk) — 1)^^ — 2xk + 2/3, to obtain 



Mi/2T y 



9^^ ^ l^sinh(Mi/2T) ] 



(53) 



so that the overall factor coincides with the one in eq. (|4^ ) for (e*). 



To estimate the goodness of the approximate result in eq. (^ we evaluated the term 
in square brackets there, taking = 5 for definitenes, and we obtained 1.163, 1.055 
and 1.000 for M\jT = 1, 3 and 10 respectively. A numerical evaluation of the exact 



correction due to the velocity, using eq. (^Tj), leads for the same factor the values 1.222, 
1.052, 1.000 respectively, so that again the accuracy is reasonable, and we see that for 
T < Ml the velocity dependent correction to e is not large. 



5 Conclusion 



The baryon asymmetry of the Universe is probably the most important manifestation 
of the existence of CP violation. In the classic scenarios for baryogenesis through the 
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out of equilibrium decay of heavy particles in the early stages of rapid expansion of 
the Universe, the asymmetries in the B (or L) violating decay rates to conjugate final 
states arise at the one loop level, involving the virtual exchange of light particles, such as 
quarks, leptons or Higgs bosons. We have studied in this paper the effects of the thermal 
background of standard particles present during the decay epoch in the evaluation of 
the CP violating asymmetries. We first reconsidered the triplet scalar decay in SU{5), 
finding that the asymmetry is reduced (contrary to an earlier result) , as could be expected 
on the basis of the Pauli exclusion principle applied to the virtual fermionic lines in the 
loop. We also included the CP violation produced by the mixing among different heavy 
states. Confronting with the T — results, the modification produced by the thermal 
effects could be as large as 50% for T ~ Mi, but diminishes with decreasing temperature, 
becoming negligible for T < Mi/ 10. 

The realization that standard model anomalous B and L violating processes are in 
equilibrium at temperatures above the electroweak phase transition one, has made the 
SU{5) scenario just of academic interest, since no net B — L asymmetry (the only one 
unaffected by anomalous processes) is generated within it. However, the same anomalous 
processes have allowed some new very attractive possibilities, including the baryogenesis 
at the electroweak scale itself (although its practical implementation faces several diffi- 
culties). The most simple and promising scenario seems to be the leptogenesis, in which 
heavy right handed neutrinos generate a lepton asymmetry in their decay, which is then 
reprocessed into a baryon asymmetry by the standard model anomalous processes. We 
showed that in these scenarios the leading thermal corrections canceled among them- 
selves, due to the opposing effects produced by the bosons and fermions involved in the 
loop. Also in the supersymmetric version of leptogenesis the thermal corrections to the 
heavy scalar neutrino decay were shown to vanish, once the thermal modification of the 
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branching ratios to the different final states are included. In view of this, we studied 
the correction due to the fact that in general the decaying particle is not at rest in the 
thermal bath. We showed that the CP asymmetry depends on the particle motion, since 
the background density distribution, and hence the thermal corrections, are now modi- 
fied in the rest frame of the decaying particle. This effect can however change the decay 
asymmetries only by at most ~ 20% with respect to the usual T = results, and hence 
these last can be safely employed. 

Acknowledgements 

This work was supported in part by CICYT under grant AEN-96/1718 and DGICYT 
under grant PB95-1077 (Spain), and by EEC under the TMR contract ERBFMRX- 
CT96-0090. 

References 

[1] K. Takahashi, Phys. Rev. D29 (1984) 632. 

[2] A. Yu. Ignatiev, V. A. Kuzmin and M. E. Shaposhnikov, JET? Lett. 30 (1979) 
688. 

[3] F. J. Botella and J. Roldan, Phys. Rev. D44 (1991) 966. 

[4] J. Liu and G. Segre, Phys. Rev. D48 (1993) 4609. 

[5] L. Covi, E. Roulet and F. Vissani, Phys. Lett B384 (1996) 169. 

[6] J. Liu and G. Segre, Phys. Rev. D49 (1994) 1342. 

[7] L. Covi and E. Roulet, Phys. Lett. B in press, [hep-ph/9611425| . 

17 



[8] A. Yildiz and P. H. Cox, Phys. Rev. D21 (1980) 906. 

[9] D.V. Nanopoulos and S. Weinberg, Phys. Rev. D20 (1979) 2484. 

[10] M. Le Bellac, Thermal Field Theory, Cambridge University Press, 1996; N. P. 
Landsman and Ch. G. van Weert, Phys. Rep. 145 (1987) 141. 



[11] R. Kobes, Phys. Rev. Lett. 67 (1991) 1384; M. E. Carrington and U. Heinz, [hep^ 
th/ 96060551 . 



[12] R. Kobes and G. W. Semenoff, Nud. Phys. B272 (1986) 329; F. Gehs, [hep= 



ph/9701410 



[13] N. S. Manton, Phys. Rev. D28 (1983) 2212; V. A. Kuzmin, V. A. Rubakov and 
M. E. Shaposhnikov, Phys. Lett. B155 (1985) 36. 

[14] M. Fukugita and T. Yanagida, Phys. Lett. B174 (1986) 45. 



[15] M. A. Luty, Phys. Rev. D45 (1992) 455; M. Pliimacher, |hep-ph/9604229| ; M. Flanz 



E. A. Paschos and U. Sarkar, Phys. Lett. B345 (1995) 248; Errata: ibidem B382 
(1996) 447, ibidem B384 (1996) 487. 

[16] B. A. Campbell, S. Davidson and K. A. Olive, Nucl. Phys. B399 (1993) 111. 

[17] H. Murayama, H. Suzuki, T. Yanagida and J. Yokoyama, Phys. Rev. Lett. 70 
(1993) 1912. 

[18] H. A. Weldon, Phys. Rev. D26 (1982) 2789. 



18 




(a) (b) 

Figure captions: 

Figure 1: Diagrams which interfere with the tree one to produce the CP violation in the 
heavy particle decay. Fig. la gives the so called vertex contribution while Fig. lb gives 
the wave function one. 
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